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Widespread multidrug resistance in Escherichia coli has necessitated the reintroduction of older antibiotics, such as nitrofuran- 
toin. However, mechanisms by which resistance to nitrofurantoin emerges in E. coli are not well elucidated. Toward this aim, we 
sequenced two nitrofurantoin-sensitive E. coli sequence types (ST540 and ST2747) and their four nitrofurantoin-resistant deriv- 
atives generated in vitro under aerobic and anaerobic growth conditions. 
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Escherichia coli is the most common cause of urinary tract infec- 
tions (UTI) in the community (1). Multidrug resistance in 
E. coli has necessitated the reintroduction of older antibiotics, like 
nitrofurantoin. Nitrofurantoin acts by undergoing reduction by 
bacterial nitroreductases to generate toxic derivatives that attach 
to ribosomal proteins (2), causing defective transcription and 
translation in bacteria (3). First described in the 1970s, nitrofu- 
rantoin resistance in E. coli involves loss-of-function mutations in 
two genes encoding nitroreductases, nfsA and nfsB (4, 5). Al- 
though these mechanisms are well known, no previous studies 
have focused beyond the nfs genes. 

We generated nitrofurantoin-resistant isolates in vitro un- 
der aerobic and anaerobic conditions from two nitrofurantoin- 
sensitive E. coli strains (multilocus sequence types [MLST] 
ST540 and ST2747, for which MICs of nitrofurantoin are 16 
and 4 /xg/ml, respectively) derived from stool samples from 
two Belgian outpatients with UTI. The strains were subjected 
to three stepwise platings on Mueller-Hinton agar supple- 
mented with increasing nitrofurantoin concentrations (0.5- to 
4-fold MIC for the parent strain). Whole-genome sequencing 
(Pacific Biosciences, Menlo Park, CA, USA) was done on the 
two parental and four nitrofurantoin-resistant strains gener- 
ated under aerobic conditions (strains ST540-A and ST2747-A, 
with nitrofurantoin MICs of 256 and 128 /j,g/ml, respectively) 
and anaerobic conditions (strains ST540-AN and ST2747-AN, 
with nitrofurantoin MICs of 64 and 32 ju-g/ml, respectively). 
Genomic DNA was isolated with the MasterPure complete 
DNA and RNA purification kit (Epicentre, Madison, WI, 
USA), according to the manufacturer's protocol. Library prep- 
aration and sequencing reactions were performed using the 
PacBio DNA template prep kit 2.0 (3 kb to 10 kb) and the 
PacBio DNA sequencing kit 2.0 with C2 chemistry. Sequence 
runs for six single-molecule real-time (SMRT) cells were per- 
formed on the PacBio RS II sequencer with a 1 X 180-minute 
movie time/SMRT cell. The SMRT Analysis portal version 2.0 



was used for filtering the reads and subreads, with default pa- 
rameters, and postfiltered data of -320 Mb on each cell/per 
strain were considered for assembly. The six genomes were 
assembled using the Hierarchical Genome Assembly Process 
(HGAP), which is available with the SMRT Analysis packages 
and accessed through the SMRT Analysis Portal version 2.0. All 
six strains were sequenced with ~28-fold coverage spanning 
-14 scaffolds for each strain. The scaffolds were ordered 
against an in silico whole-genome map of E. coli ATCC 8739 
(GenBank accession no. CP000946) using the MapSolver soft- 
ware (Whole Genome Mapping; OpGen, Gaithersburg, MD) to 
derive a complete gapless chromosome (B. B. Xavier, J. Sabi- 
rova, P. Moons, J. -P. Hernalsteens, H. De Greve, H. Goossens, 
and S. Malhotra-Kumar, submitted for publication) . The chro- 
mosome sizes of the ST540 strains are 4,758,628 bp, 4,807965 
bp, and 4,875,674 bp. The ST2747 strains have chromosome 
sizes of 5,054,424 bp, 4,998,910 bp, and 5,090,326 bp, with a 
G-l-C content of -51%. We used the Rapid Annotations using 
Subsystems Technology (RAST) server (6) and NCBI Prokary- 
otic Genome Automatic Annotation Pipeline (PGAAP) (7) for 
function-based annotation for all strains. The genomes of 
ST540 and ST2747 contain 4,587 and 4,859 protein-coding 
genes and 114 and 119 RNA genes, respectively, of which 89 
encode tRNAs and 8 encode rRNAs. In addition, our genome 
annotation confirmed the duplication of the mobile elements 
and other single-nucleotide polymorphisms (SNPs) in ST540 
and ST2747, which were identified using an inbuilt tool in the 
CLC Genomics Workbench 6.5.1 (CLC, Inc., Aarhus, Den- 
mark). 

Nucleotide sequence accession numbers. The genome se- 
quences for all sbc strains have been deposited at DDBJ/EMBL/ 
GenBank under the accession no. CP007265, CP007390, 
CP007391, CP007392, CP007393, and CP007394. The versions 
described in this paper are the first versions. 
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